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Abstract 

We report the observation of the decay (jxpK^ and find evidence for cfxpK^. These 

results are based on a 414 fb~^ data sample that contains 449 x 10^ BB pairs, collected with 
the Belle detector at the KEKB asymmetric-energy e"'"e~ (3.5 on 8 GcV) collider operating at 
the T(45') resonance. This is the first observation of a 6 ^ s'ss'ss transition. The branching 
fractions for these decay modes are measured to be B{B^ 4>(j)K^) = (3.2lo;5 it 0.3) x 10~^ 
and ^ = (2.3tJ:? ± 0-2) x IQ-^ for 00 invariant mass below 2.85 GeV/c^. The 

corresponding partial rate asymmetry for the charged B mode is measured to be Acp{B^ — 
(jxpK^) = COl^Q^g lb 0.02. Results for other related charmonium decay modes are also reported. 
We also search for CP asymmetry using the 00 candidates within the rjc mass region. The value 
obtained is Acp{B^ ^ <p<pK^,M^^ w M^J = 0.15lo;}7 ± 0.02, which is consistent with no 
asymmetry. 

PACS numbers: 13.25.Ft, 13.25.Hw, 14.40.Nd 
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Evidence of charmless B (jxpK decays has been reported by the Belle collaboration 
using 85 X 10^ BB pairs [l^. In the standard model (SM), this decay channel requires the 
creation of an additional final ss quark pair in a 6 — >■ sss process, such as i? — (pK. Although 
the statistical error of the previous B — >■ (jxpK measurement was large, the branching fraction 
was found to be around 1/3 of that oi B ^ (pK, which provided useful information for 
understanding quark fragmentation in B decay. Moreover, it indicated that with sufficient 
statistics the decay B (pipK could be used to search for a possible non-SM CP-violating 
phase in the h ^ s transition [2j. Direct CP violation could be enhanced as high as the 
40% level if there is sizable interference between the transitions due to non-SM physics and 
the decays via the rjc resonance. 

In this paper, we report the observation of charmless and charmful B (jxpK decays 
based on a 414 fb~^ data sample that contains 449 xlO^ BB pairs, collected with the Belle 
detector at the KEKB asymmetric-energy e+e" (3.5 on 8 GeV) colhder g operating at the 
T(45') resonance. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-fiight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux- 
return located outside of the coil is instrumented to detect K° mesons and to identify muons 
(KLM). The detector is described in detail elsewhere 0]. Two inner detector configurations 
were used. A 2.0 cm radius beampipe and a 3-layer silicon vertex detector (SVDl) were used 
for the first sample of 152 xlO^ BB pairs, while a 1.5 cm radius beampipe, a 4-layer silicon 
detector (SVD2_) and a small-cell inner drift chamber were used to record the remaining 297 
X 10*^55 pairs [|. 

Charged kaons are required to have impact parameters within ±2 cm of the interaction 
point (IP) along the z-axis (anti-parallel to the positron direction) and within 0.2 cm in the 
transverse plane. Each track is identified as a kaon or a pion according to a K/n likelihood 
ratio, 7l{K/7i) = Ck/{Ck + C-,,), where CkI is the likelihood of kaons/pions derived from 
the responses of TOF and ACC systems and the energy loss measurements from the CDC. 
The likelihood ratio is required to exceed 0.6 for kaon candidates; within the momentum 
range of interest, this requirement is 88% efficient for kaons and has a misidentification rate 
for pions of 8.5%. Neutral kaons are reconstructed via the decay A'^ — > 7r"'"7r~ and have an 
invariant mass 0.482 GeV/c^ < M^^+t^- < 0.514 GeV/c^ (±4(T mass resolution). The 7r"'"7r~ 
vertex is required to be displaced from the IP and the fiight direction must be consistent with 
a Kg that originated from the IP. The required displacement increases with the momentum 
of the Kg candidate. 

B meson candidates are reconstructed in the five kaon final state. Two kinematic variables 
are used to distinguish signal candidates from backgrounds: the beam-energy constrained 

mass Mbc = \/-Ebeam " l-PrcconP and the energy difference A^; = E.^con - -^beam, where -Ebeam 
is the beam energy, and -Erecon and -Precon are the reconstructed energy and momentum of 
the signal candidate in the T(4S') rest frame. The resolution of Mbc is approximately 2.8 
MeV/c^, dominated by the beam energy spread, while the AE resolution is around 10 MeV. 
Candidates with five kaons within the region |Aii^| < 0.2 GeV and 5.2 GeV/c^ < Mbc are 
selected for further consideration. The signal region is defined as 5.27 GeV/c^ < Mbc < 5.29 
GeV/c2 and \AE\ < 0.05 GeV. 

The dominant backgrounds are e~^e~ ^ qq [q = u, d, s, c) continuum events. Event 
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topology and B flavor tagging are used to distinguish the jet-like continuum events and the 
spherically distributed BB events. Seven event-shape variables are combined into a single 
Fisher discriminant j|6ji|. The Fisher variables include the angle between the thrust axis of 
the B candidate and the thrust axis of the rest of the event (cos 6*7-), five modified Fox- 
Wolfram moments Q, and a measure of the momentum transverse to the event thrust axis 
{S±) Two other variables that are uncorrelated with the Fisher discriminant and help 
to distinguish signal from the continuum are cos 6b, where is the angle between the B 
flight direction and the beam direction in the e~^e~ center of mass frame, and Az, the z 
vertex difference between the signal B candidate and its accompanying B. We form signal 
and background probability density functions (PDFs) for the Fisher discriminant, cos9b 
and Az using the signal Monte Carlo (MC) events and sideband data (5.2 GeV/c^ < Mbc < 
5.26 GeV/c^), respectively. The products of the PDFs for these variables give signal and 
background likelihoods Cs and Cbg for each candidate, allowing a selection to be apphed 
to the likelihood ratio 71 = Cs/ {Cs + Cbg)- 

Additional background discrimination is provided by the quality of the B flavor tagging 
of the accompanying B meson. The standard Belle flavor tagging package P] gives two out- 
puts: a discrete variable indicating the flavor of the tagging B and dilution factor r, which 
ranges from zero for no flavor information to unity for unambiguous flavor assignment. The 
continuum background is reduced by applying a selection requirement of TZ for events in 
each r region according to the figure of merit defined as Ns / V^s + ^bg, where Ns denotes 
the expected (jxpK signal yields based on MC simulation and the branching fraction reported 
in our previous measurements, and Nbg denotes the expected qq yields from sideband data. 
This requirement removes (61-81)% of the continuum background while retaining (80-92)% 
of the signal, and depends on the decay channel {(pipK or (jxpK^) and the SVD configu- 
ration during the measurement (SVDl or SVD2). Backgrounds from other B decays are 
investigated using a large MC sample and are found to be negligible after the TZ requirement. 

The signal yields are extracted by applying an unbinned extended maximum likelihood 
(ML) fit to the events with Mbc > 5.2 GeV/c^ and \ AE\ < 0.2 GeV. For the ([xpK^ mode, we 
simultaneously obtain the yield and the partial rate asymmetry Acp using the likelihood, 
defined as: 



C = exp[-{Ns + Nb)] n(E [1 - ■ ^'cp]A^.^(Mbc, AE)), (1) 

i 3 ^ 

where i is the identifier of the i-th event, j indicates signal or background, P(Mbc, AE) is 
the two-dimensional PDF of Mbc and AE", and q indicates the B meson flavor, +1 for B^ 
and —1 for 5", respectively. For neutral B events, ^[1 — gj -^cp] in Eq- © is replaced b y 1 . 
The Mbc PDFs are modeled by a Gaussian function for signals and an ARGUS function [10| 
for the continuum, while a Gaussian is used to describe the signal AE and a 2"'^ order 
Chebyshev polynomial is used for the background AE. The parameters of the PDFs are 
determined using high-statistics MC samples and sideband data for signal and background 
shapes. The signal PDFs are calibrated by comparing the Mbc and AE distributions of the 
B^ D^{K^ 71^ 71^ 11^)71^ events with the MC expectation. 

We search for charmless B (fxpK decays by requiring the (jxj) invariant mass (M^^) to 
be less than 2.85 GeV/c^ the region below the charm threshold. Candidate mesons are 
identified by requiring the invariant masses of K^K~ pairs {Mk+k-) to be in the range 
1.0 GeV/c^ to 1.04 GeV/c^ (±4.6cr). Figure^shows the Mbc and AE projections with the 
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fit curves superimposed. Clear signals appear in both and modes with signal yields 
of 37.Q11qI and T.Slil'^, respectively. Although K^K~ candidates are required to lie in the 
(j) mass region, non-0 backgrounds may also contribute. Figure |2ta) shows the Mk+k- vs. 
Mk+k~ distributions for {K^ K~ K~)K^ candidates in the signal region, where the two 
K^K^ pairs are required to have invariant masses less than 1.2 GeV/c^. Events in the 
two bands are used to estimate the B^ (pK^K'K^ contribution. Figure |2fb) shows 
B signal yields as a function of the K^K~ invariant mass after requiring the other 
K^K~ pair to have a mass in the mass region. The B signal yields are fitted with a 
threshold function in the region 0.98 GeV/c^ < Mk+k- < 1-2 GeV/c^, excluding mass 
region (1.0 GeV/c^ < Mk+k- < 1-04 GeV/c^). The amount of the non-0 contribution 
is estimated by interpolating the B yields in the sideband region to the mass region, 
which is 4.4lo;7 events. Since events in the two bands contain both true mesons and 
non-resonant K^K~ pairs, the area underneath the mass region in Fig. I2tb) includes the 
(j)K^ K~ K"^ contribution but counts the non-resonant bK component twice. Therefore, we 
estimate the non-resonant B — > bK contribution using the B signal yield in the upper right 
corner of the dashed region in Fig. |2^a). We assume a phase-space distribution in 4-kaon 
mass. We obtain 1.3 ± 0.4 non-resonant events in the (fxpK^ sample. After subtracting this 
contribution of 1.3 ± 0.4 events, the non-00ii' fraction is calculated to be (7 ± 4)%. The 
same procedure is applied to the (jxpK^ sample; here we obtain a fraction of (7 ± 9)%. 

Table m summarizes the (jxpK results after subtracting the non-00ii' contribution. Signal 
efficiencies are obtained by generating (fxpK MC events, where the same M^^ < 2.85 GeV/c^ 
requirement is applied. Systematic uncertainties in the fit are obtained by performing fits 
in which the signal peak positions and resolutions of the signal PDFs are successively varied 
by ±lcT. The quadratic sum of each deviation from the central value of the fit gives the 
systematic uncertainty of the fit. For each systematic check, the statistical significance 

is taken as — 2 ln(£feeddown/ -Cmax) , where £feeddown and £max are the likelihoods at the 
expected non-00ii' yields and the best fit, respectively. We regard the smallest value as our 
significance including the systematic uncertainty. The number of B^B^ and B^B^ pairs are 
assumed to be equal. 



TABLE I: Signal yields, efficiencies including secondary branching fractions, significances, branch- 
ing fractions for B (j)(f)K and related charmonium decays. 

Mode 

B± ^ ^<^i^± (M^0 < 2.85 GeV/c2) 
^ 00^0 (M^^ < 2.85 GeV/c2) 
ricK^, Vc — > (t>4> 

B± ^ J/^pK^, J/tp (j)K+K- 
J/^K^, J/i}} -> 2{K+K-) 



Yields 


efBciencies(%) 


E 


i3(10-6) 


34.21^;^ 


2.41 


9.5 


3.2t°j±0.3 


7 q+3.0 
' -"^-2.4 


0.69 


4.7 


2.3tJ;?±0.2 




2.72 


7.2 


2.4t°j±0.2 


76.8tl^i 


4.85 


9.4 


3.5 ±0.6 ±0.3 


104.6ti?;3 


9.93 


10.2 


2.4+°:^ ±0.2 




4.67 


8.5 


1.2 ±0.3 ±0.1 




9.41 


9.7 


0.97l^;i^±0.1 



The performance of the IZ requirement is studied by checking the data-MC efficiency 
ratio using the B^ D^{-^ 71^71^ 71^)71^ sample. The obtained error is 2.7-2.8%. The 
systematic errors on the charged track reconstruction are estimated to be around 1% per 
track using partially reconstructed D* events. Therefore, the tracking systematic error is 
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FIG. 1: Distributions of Mbc and AE for the decay modes 4>(pK^ (a,b) and 5° (pcpK^ 

(c,d), with (jxp invariant mass less then 2.85 GeV/c^. The open histograms represent the data, 
the soUd blue curves show the result of the fit, the dash-dotted red lines represent the signal 
contributions and the dashed green curves show the continuum background contributions. 



5% (5 tracks) for the (jxpK^^ mode and 4% for the (jxpK^ mode. The kaon identification 
efficiency is studied using samples of inclusive D*^ D^ti^^D^ K~ti^ decays. The 
K'g reconstruction is verified by comparing the ratio of — >■ K^tt^ and K^ti^tt^ 
yields. The resulting detection systematic error is 4.9%. The uncertainty in the number 
of BB events is 1%. The final systematic error is obtained by summing all correlated errors 
linearly and then quadratically summing the uncorrelated errors. 

After subtracting the non-^^i^' contribution, the branching fractions for charmless B — >■ 
(p<pK decays are B{B^ (f)(t)K^) = (3.2l|]:^±0.3) x IQ-^ with 9.5a significance and B{B^ 
(j)(j)K^) = (2.3^0 7 ± 0.2) X 10~^ with 4.7(T significance. The measured charge asymmetry for 
5± ^ 4>(^K^ decay is O.Oll[!;J^ ± 0.02. The first error is statistical and the second is 
systematical. 

It is of interest to search for possible 00 resonances above charm threshold. Figure Efa) 
shows the B signal yields divided by the bin size as a function of M,^^ after releasing the 
M^fP < 2.85 GeV/c^ requirement. There is no specific enhancement in high 00 mass region 
except for the r]c peak near 3 GeV/c^. Reference Q suggests a large CP asymmetry, that 
arises from the interference between i?^ — > (fxpK^ and B^ Vc{~>' 00)-?^^ decays. The 
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FIG. 2: (a) The distribution of M^+r- vs. M^+k- for the K+K-R+K-R^ candidates in the 
Mbc — Ai? signal box with Mr+r- < 1.2 GeV/c^. The two R^R^ bands indicate the (j) mass 
region (1.0 GeV/c^ < M^+j^- < 1.04 GeV/c^). The horizontal and vertical dashed lines are 
located at M^+^- = 1.05 GeV/c^. The rectangle on the upper right is the sideband region; 
events in this region are used to estimate the non-resonant B 5R contribution, (b) B signal 
yields as a function of the Mj^+j^- after requiring the other R^R^ pair has a mass in the (p mass 
region. The threshold function is used to be the fit curve and events with 1.0 GeV/c^ < Mr+r- < 
1.04 GeV/c^ are excluded when the fit is performed. 

events with (jxj) invariant mass within ±40 MeV/c^ of the nominal rjc mass are selected 
to investigate this asymmetry. The measured CP asymmetry is 0.15!tai7 ± 0.02, which is 
consistent with no asymmetry. 

We study possible charmonium states by checking B yields with M^k between 2.8 GeV/c^ 
and 3.2 GeV/c^. Since rj^ and J/ip mesons may decay to cj)K^K~ and 2{K^K~) pairs, a 
mass scan is performed with and without the requirement that the K^K~ pair lie in the 
mass region. As shown in Fig. El clear rjc and J/ip resonances are visible in the (j)K~^K~ 
and AK samples while only an rjc peak appears in the 00 mode. 

We obtain the signal yields for f]cK^ and B^ J/ipK^ by performing binned 

histogram fits to Figs. intb),infc) andEfd). The J/ip signal PDF is modeled with a Gaussian 
function while the rjc PDF is described by a Breit-Wigner function convolved with a Gaussian 
resolution function, which has the same Gaussian width as the J/ip PDF. Since sizable signals 
are observed in the AK mode, the parameters are determined using the 4K sample and the 
same signal PDFs are then applied to the (j)K^K~ and 00 samples. The obtained Gaussian 
width is measured to be 4.0tJ;8 MeV/c^ by performing a fit with a 2"°' order Chebyshev 
polynomial as the non-resonant PDF after excluding events in the rjc mass region (2.94 
GeV/c^ < M^^K < 3.02 GeV/c^). A fit is performed to the full range with the Gaussian 
width fixed. The signal yields are summarized in TablelH The peak positions obtained for the 
T]c and J/i) are 2.979 ±0.002 GeV/c^ and 3.094 ±0.001 GeV/c^, respectively, consistent with 
the nominal r]c and J /ip masses. The rjc Breit-Wigner width is measured to be 25.2t6;o i 0-3 
MeV/c^, where the central value is consistent with the world average and the second error 
is due to the uncertainty in the mass resolution. 

For the (j)K~^K~ and 00 modes, the non-0 contribution is determined from the B signal 
yields for events with one K^K~ pair in the sideband region (1.05 GeV/c^ < Mr+r- < 
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FIG. 3: signal yields as a function of (a,b) M^^, (c) M^^+j^- and (d) M4x- In (a) we use 
different bin sizes for Mfp^ less than 3 GeV/c^ and greater than 3 GeV/c^. The subset with M^^ 
from 2.8 GeV/c^ to 3.2 GeV/c^ is shown in (b). 

1.09 GeV/c^) and the 4if and ipK^K^ masses are in the charmonium resonance region, 
respectively. We find 3. 01^;° events in the rfc — > 4>4> mode, 6.4^4;! events in the 77^ (j)K^K~ 
mode, and 3.4110 in the J /ip ^ (j)K~^K~ mode. After subtracting the feed-down yields, we 
obtain the results listed in Table HJ 

Signal efficiencies are determined using signal MC and the detection systematic uncertain- 
ties are similar to what was described in the charmless (fxpK part. Fit systematic uncertain- 
ties are estimated by successively varying the peak positions and resolutions of the Mbc — AE' 
signal PDFs as well as the convolution Gaussian width in the fit. The quadratic sum of each 
deviation gives the fit systematic errors. Since the sub-decay branching fractions of rjc and 
J/ ip mesons to AK, (pKK and (jxp final states are not precisely known, we provide the product 
of branching fractions for various decays in Table HJ Using the known branching fractions of 
B{B^ 7]cK^) = (9.1 ± 1.3) X 10"^ and B{B^ J/tpK^) = (1.008 ± 0.035) x 10"=^ [l^, 
the subdecay branching fractions are calculated and listed in Table El 

In summary, we have observed the charmless decay — > (fxpK^ and evidence of 
B^ (j)(f)K^. We also report the CP asymmetry of the charged decay and measurements 
of other closely related charmonium decays. The results are consistent with the previous 
measurements, but have considerably improved precision due to the increase in statistics. 
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TABLE II: The measured branching fractions of secondary charmonium decays and the world 
averages [l^ . 



Decay mode 



B (measured) 



B (PDG 2006 value) 



r]c (t>K+K- 
T], ^ 2{K+K-) 
Jji) (l)K+K- 
J/ij 2{K+K-) 




(2.7 ±0.9) X 10-3 
(2.9 ± 1.4) X 10^3 
(1.5 ±0.7) X 10-3 
(1.83 ±0.24) X 10-3 
(7.8 ± 1.4) X 10-^ 
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